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We study the joint spectral properties of photon pairs generated by spontaneous parametric 
down-conversion in a one-dimensional nonlinear photonic crystal in a collinear, degenerate, type- 
II geometry. We show that the photonic crystal properties may be exploited to compensate for 
material dispersion and obtain photon pairs that are nearly factorable, in principle, for arbitrary 
materials and spectral regions, limited by the ability to fabricate the nonlinear crystal with the 
required periodic variation in the refractive indices for the ordinary and extraordinary waves. 
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I. INTRODUCTION 

Two-photon states with specific continuous variable- 
entanglement properties are required for a number of 
quantum-information-processing applications. In partic- 
ular, pure-state single photons, crucial for applications 
relying on interference between independently prepared 
single photons, such as quantum computing with lin- 
ear optics [1], entanglement swapping [2J, and tclcpor- 
tation [3], can be heralded only from factorable photon 
pairs, where no correlations exist between the constituent 
single photons [J. Let us note that factorability may be 
imposed by postselecting an appropriate subensemble of 
photon pairs through spectral and/or spatial filtering; 
this can be done only at the cost of a typically dras- 
tic reduction in count rate and source brightness. Scal- 
able quantum-information processing requires the elimi- 
nation of this type of postselection. It has furthermore 
been shown that, for quantum- information applications 
involving multiple-pair generation resulting from high- 
gain parametric down-conversion, spectral or spatial fil- 
tering is likewise ineffectual. [S]. In addition, the coex- 
istence of factorability and a large generated bandwidth 
is important for some applications, such as the genera- 
tion of ultrashort Fourier-transform-limited heralded sin- 
gle photons [B] . The process of spontaneous parametric 
down-conversion represents a well-established method for 
the generation of photon pairs, leading to unparalleled 
flexibility in the resulting continuous-variable entangle- 
ment properties. Indeed, the nonlinear crystal dispersion 
characteristics in conjunction with temporal and spatial 
structure in the pump beam may be exploited to engineer 
the type and degree of correlations present in spectral and 
transverse wave-vector degrees of freedom at the source, 
thus eliminating the need for photon pair filtering. 

The use of a broadband pump is essential in the con- 
text of the generation of photon pairs with spectrally 
engineered properties; indeed, a monochromatic pump 
can access only a one-dimensional subspace of signal- 
idler frequency space, which precludes certain states of 
interest. Previous work has established the central role 
that is played by group velocity mismatch terms between 



the three fields involved in parametric down-conversion 
(pump, signal, and idler) in the determination of the re- 
sulting photon pair properties. Thus, in Rcf. [7], it was 
shown that if the pump pulse propagates at a higher 
group velocity than one of the generated photons but 
lower than the conjugate generated photon, then it be- 
comes possible to emit factorable photon pairs. An im- 
portant limitation of this technique is that the group 
velocity condition which must be satisfied occurs only 
for specific materials, at specific spectral ranges [H [1], 
typically in the infrared where single-photon detection is 
unfortunately not well developed. A number of tech- 
niques have been proposed and in some cases imple- 
mented which enable effective control over the photon 
pair entanglement properties in the spectral domain un- 
constrained by material dispersion. In one such tech- 
nique, the effective group velocities experienced by the 
three fields involved are modified by a pair of diffrac- 
tion gratings placed before and after the nonlinear crys- 
tal so as to generate photon pairs with tunable spec- 
tral characteristics (TUJ |TT]. Alternatively, the spectral 
content of two photon states may be restricted to the 
modes supported by a nonlinear cavity which leads, for 
a short, high-finesse cavity, to factorable, narrowband 
photon pair generation |12j . Likewise, on the one hand a 
transversely pumped source where signal and idler pho- 
tons are emitted in a counterpropagating waveguided ge- 
ometry |13j and on the other hand noncollinear paramet- 
ric down-conversion (PDC) where a specific relationship 
between transverse and longitudinal phase matching is 
attained [H] can both lead to states with spectrally en- 
gineered properties. Another route is the use of crys- 
tal superlattices, where the dispersion in short crystal 
segments is compensated by birrefringent compensators, 
permitting two-photon states with a wide range of pos- 
sible spectral entanglement characteristics |T5] . 

Previous theoretical work has explored the use of non 
linear photonic crystals in the context of the process 
of spontaneous parametric down-conversion. In partic- 
ular, Refs. [inillT] study the potential of semiconductor- 
based nonlinear one-dimensional photonic crystals to 
yield phase-matching properties appropriate for the gen- 
eration of polarization-entangled photon pairs. Likewise, 
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it has been shown that even a short one-dimensional 
photonic crystal is capable of generating a considerable 
photon pair flux due to field localization in such struc- 
tures [T^ [ini dO] . In this paper we study the potential 
of exploiting the properties of one-dimensional nonlinear 
photonic crystals in order to generate photon pairs with 
engineered spectral entanglement properties. We con- 
centrate on type-II, frequency degenerate, coUinear PDC. 
Collinear operation is important because it permits PDC 
in a waveguided geometry, which leads to larger genera- 
tion rates as well as to effective control over the transverse 
spatial structure of the emitted light, crucial for efficient 
single-mode fiber coupling and for spatial mode match- 
ing in interference experiments |21j . Type-II operation 
for which the signal and idler photons are orthogonally 
polarized is important because it enables efficient photon 
pair splitting. We show that a weak index of refraction 
modulation present in an otherwise standard birefringent 
nonlinear crystal can be exploited to attain the group 
velocity conditions, in addition to basic phase matching, 
required for factorable photon pair generation. 



II. ONE-DIMENSIONAL NONLINEAR 
PHOTONIC CRYSTALS 
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FIG. 1: (color online) Schematic of a one-dimensional, non- 
linear photonic crystal, (a) Second order nonlinearity. (b) 
Refractive index, shown for one of the polarizations, (c) Rep- 
resentation of the periodic material, with h = h. — a. 



In this paper we analyze the generation of photon pairs 
by the process of parametric down-conversion in a x(2) 
material characterized by a spatial periodicity in its lin- 
ear optical properties. In particular, we study PDC pro- 
duced by a nonlinear photonic crystal (NLPC) based on 
a standard bulk nonlinear crystal with uniaxial birefrin- 
gence which has been modified from its natural state by 
an appropriate physical mechanism so as to yield a pe- 
riodic variation of the index of refraction for both po- 
larizations, while maintaining the nonlinearity constant 
throughout the crystal thickness. Concretely, within one 
period of the resulting Bragg grating, we assume that 
the index of refraction for the ordinary (/i = o) and the 
extraordinary ray (/i = e) are given by 

= / "'^^H' 0<^<«' (1) 
' ^ \ n^2{<^), a<z <K. > 

This unit cell is assumed to be replicated throughout 
the crystal length. Here, A represents the Bragg period 
and a/ A is the duty cycle. We refer to such a material 
as a one-dimensional nonlinear photonic crystal. Fig- 
ure [T] shows a crystal schematic where A indicates zones 
characterized by indices of refraction Uqi and n^i and B 
indicates zones characterized by indices of refraction no2 
and 11^2- We assume that in zones A the crystal has been 
left in its natural state, while in zones B it has under- 
gone index change. For simplicity, we assume that zones 
A and B are of equal length (leading to an a/A = 0.5 
duty cycle) and we likewise assume that the permittiv- 
ity contrast a, equivalent to the index square contrast, is 
frequency independent. 



^ ^ 2[e^i(a;) - £p2(c^)] ^ ~ "mzH^ 

where 

n^{ojf = Ki(co)2 ^ „^2H']/2. (3) 

with /i = e, o. Throughout this paper we will assume 
that propagation of the three fields involved is normal to 
each of the interfaces between zones A and B. Under 
these conditions, such a material can exhibit, for each of 
the polarizations, a so-called photonic band gap centered 
at each of the Bragg wavelengths 

y^Bragg,ra ^ 2n^A/m, (4) 

where m = 1,2,3,..., and where /i = e,o. Within each 
band gap, for a sufficient crystal thickness, light is effi- 
ciently reflected, while for frequencies outside the band 
gap, light propagates in the form of so-called Bloch waves 
which can be written as 

E{z, t) = Ek{z) exp{i[K{uj)z ~ cut]}. (5) 

Here, Ej({z) is the Bloch envelope, which exhibits a 
spatial periodicity matching that of the material, i.e., 
Ek{z -I- A) = Ek{z), while K{u!) represents the Bloch 
wave number. Note that, for a continuous material with- 
out Bragg grating, the Bloch envelope reduces to a con- 
stant, and therefore the Bloch wave reduces to a plane 
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wave. Following a coupled-mode analysis where the spa- 
tial periodicity in the permittivity is assumed to be well 
represented by a weak perturbation to the material per- 
mittivity, it can be shown that the Bloch wave number 
in the vicinity oi K = rmr/A can be expressed as [22] 



4™)(c.) = nm/A ± V^(A/3^"V2)2-|4"^P 



(6) 



Here, k^™"* represents the coupling coefficient between the 
forward and backward waves 



v^™^ = »[1 - cos(m7r)]afc^(a;)/(47rm). 



(7) 



and A/?^"*-* represents the Bragg phase mismatch between 
these two waves, 



(8) 
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FIG. 2: (Color online) Plots as a function of frequency of 
(a) the Bragg wave number Kiu)), (b) the group velocity 
IjK'iiJ)^ (c) K'(lS)^ and (d) group velocity dispersion co- 
efficient K"{lj); for comparison, tfie dotted horizontaf fine 
indicates the magnitude of the GVD coefficient for SF-10, a 
particufarfy dispersive giass. 



where fc^ — n^u/c characterizes the underlying material 
dispersion; note that Eq.Q follows from the condition 
A/3(") = 0. A photonic band gap appears when k^™-* ^ 
and its effects become appreciable for A/?^™'' w 0. The 
band gap boundaries Umin and LUmax are obtained by 
solving for the frequencies such that the argument of 
the square root in Eq.([6]) vanishes. For a forward- 
propagating wave, the minus sign in front of the square 
root in Eq.([6| should be used for w < tOmin while the 
plus sign should be used for to > ujmax- Let us note that 
for a 0.5 duty cycle, Eq.Q tells us that, even order band 
gaps (m = 2,4,6...) are suppressed. 

In what follows, we will concentrate on first order 
(to = 1) band gaps, though the analysis below could 
be employed for any order. In particular, for practical 
implementations it may be challenging to fabricate the 
required Bragg gratings with periods in the region of hun- 
dreds of nanometers compatible with modified dispersive 
characteristics in the visible; alternatively, it is possible 
to use longer periods, together with higher-order band 
gaps. The dispersive properties of NPLCs are charac- 
terized by the function K^to), where from this point we 
omit the order superscript. For propagating waves at 
frequencies in the vicinity of the band gap boundaries, 
K(uj) can differ substantially from the underlying mate- 
rial dispersion relation k{iv). In particular, group veloci- 
ties can be strongly reduced, while group velocity disper- 
sion terms (as well as higher-order frequency derivatives 
of the wavenumber) can increase markedly, from their 
values associated with material-only dispersion. 

In order to illustrate these effects, let us consider a 
specific example of a one-dimensional NLPC, based on 
a /9-barium-borate (BBO) crystal. We assume that the 
Bragg period is given by A = 279.1 nm with a a/ A = 0.5 
duty cycle, the angle subtended by the pump and the 
optic axis (to be referred to as the propagation angle) is 



41.8°, and the permittivity contrast is a = 0.027. Fig- 
ure [2] shows the resulting dispersive properties. Figure [2] 
(a) shows, for each of the two polarizations, a plot of the 
Bloch wave number K(uj), showing clearly a band gap 
for each of the polarizations, within which K{lo) becomes 
complex- valued. Figure |2] (b) shows a plot of the group 
velocities for each of the two polarizations, exhibiting a 
marked reduction near the band gap boundaries. Fig- 
ure [2] (c) shows the first frequency derivative of K{u}), or 
the reciprocal group velocity. Finally, Figure [2] (d) shows 
a plot of the group velocity dispersion (GVD) term for 
each of the polarizations. It is apparent from the plot 
that GVD can be greatly enhanced near the band gap 
boundaries, and likewise that it is possible to obtain both 
positive and negative GVD. 



III. PDC IN ONE-DIMENSIONAL PHOTONIC 
NONLINEAR CRYSTALS 

Following a standard perturbative approach, the quan- 
tum state describing photon pairs produced by paramet- 
ric down-conversion in the spontaneous limit may be ex- 
pressed as 



m)) 



1 + - / dt'H{t') 
in Jo 



I vac), 



(9) 



where |vac) denotes the vacuum and H is the interaction 
Hamiltonian 



H{t)^eol dV d{r)E'^\r,t)Elr\f,t)E\ ^(f,i)+H.C. 

(10) 
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Here, V is the illuminated volume in the nonlinear 
medium, d(r) is the second-order nonlinearity, and 
Ef^{f,t) {fi — p,s,i) represents the electric field oper- 
ators associated with each of the interacting fields. In 
a nonlinear photonic crystal, each of these fields is de- 
scribed by a Bloch wave. Thus, if we assume that the 
pump field is classical, it can be expressed as 



E^+'> (f, i) -> du Up (llj) Ekj,{z, tj)exp {i [Kp{ijj)z - ujt]} 



(11) 

in terms of the Bloch wavenumber Kp{Lo), Bloch enve- 
lope EKp{z,uj), and the spectral amplitude ap{Lu). It 
is convenient to express the Bloch envelope as a Fourier 
series. 



EK^{z,Lo)=Y,£A^y'^'' 
I 



(12) 



in terms of the spatial harmonics Gi — 27r//A. The signal 
and idler fields are quantized; their positive-frequency 
part can be expressed as 



El+\r,t) = i j duoY,£^.i{^y^{^)a^{K^{uj) + Gi) 



xexp {i {[K^,{uj) + Gi]z - uot}) 



(13) 



in terms of the Bloch wave number K^{u)), envelope 
Ek^{z,uj) and the Bloch envelope Fourier series co- 

eSicients e^;. Here, £^(a)) = yJhujK'^iuj) /[2efj,{uj)S], 

where K'^(uj) is the first frequency derivative of K^, 
(XpK{ijj) is the annihilation operator for the signal(s) or 
idler(i) modes, e^(a;) is the permittivity in the nonlinear 
medium, and S is the transverse beam area. 

It can be shown that the resulting two-photon com- 
ponent of the state may be written for specific direc- 
tions of propagation (throughout this paper assumed to 
be coUinear with the pump beam) as 



1^) - E 



Imn 



duJsduJifi„in{uJs,UJi) 



xaliK{Lj,) + G^)al{K{u;,) + G„)|vac) (14) 

where the joint spectral amplitude fimni^sT^i) may be 
factored as /imn (t^s,Wi) = ap(ws-|-Wj)0im„(ws, w^). Here, 
ap{LOs + oji) represents the pump spectral envelope func- 
tion, while (j)imn{^s,^i) is a function which describes 
the phase-matching properties of the nonlinear photonic 
crystal and can be expressed as 



0im«(Ws,Wj) = epl{uJs+UJi)£l^{uJs)£*in{^t) 

x4(a;s)4(a;i)sinc {LAKimn/2) 

X exp {iLAKijnri/'2) , (15) 



in terms of the frequency-dependent phase mismatch 
AKimn adjusted by the momentum contribution due to 
the photonic crystal structure. 



AK + 2Tr{l - m~ n)/A 



(16) 



where AK — Kp — Kg — Ki. The term proportional to 
27r/A is similar to that which appears for quasi-phase- 
matched interactions and has the effect of shifting the 
spectral range where phase-matching is attained. Thus, 
the phase-matching contributions for different values of 
I — m — n, if they exist (i.e., if Akimn = 0), will tend 
to be spectrally distinct from each other. Note, however, 
that in general terms, for photonic crystal periods in the 
hundreds of nanometers, contributions with l~m — n^O 
result in a term proportional to 27r/A which will tend to 
be larger than the wave numbers for each of the three 
interacting fields, and will therefore also be larger than 
AK; this makes it difficult to achieve AKi,nn — for 
I — m — n ^ 0. 
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FIG. 3: (Color online) Normalized modulus of Fourier co- 
efficients for (a) an extraordinary wave at 425 nm, (b) an 
extraordinary wave at 850 nm, and (c) an ordinary wave at 
850 nm. 

Let us now consider the Bloch wave characteristics 
at 425 nm and 850 nm (which could represent pump 
and PDC wavelengths) in a material with the parame- 
ters specified at the end of the previous section. Fig- 
ure [3] shows the modulus of the leading Fourier series 



coefficients [see Eqs.(12l and (13l], calculated from the 
Bloch envelope in turn determined by the eigenvectors of 
the translation matrix which characterizes the periodic 
material [52]. Figure [Sj a) shows the Fourier coefficients 
for an extraordinary wave at 425 nm. Fig. |3|b) shows 
the Fourier coefficients for an extraordinary wave at 850 
nm, and Fig. ^c) shows the Fourier coefficients for an 
ordinary wave at 850 nm. Let us note that a wave at 
425nm behaves essentially like a plane wave with negligi- 
ble m 7^ terms. This is due to the fact for this specific 
material there is no band gap in the vicinity of 425 nm; 
in particular, the second order band gap is suppressed by 
employing a a/ A = 0.5 duty cycle. This specific material 
exhibits a band gap centered at 904.9 nm for the ordi- 
nary wave and at 876.3 nm for the extraordinary wave. 
A wave of either polarization at 850 nm is sufficiently 
near to the corresponding band gap so that the resulting 
dispersion relation is significantly modified (see Fig. |2|, 
yet as observed in Figs.l3jb) and (c small m = 1 (i.e., re- 
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fleeted wave) contribution and where other contributions 
are neghgible. 

In what follows, we will study the generation of PDC 
where the emission frequency is in sufficient proximity 
to the band gap boundaries so that dispersion is strongly 
modified, though sufficiently removed so that the fraction 
of light appearing in modes other than the fundamen- 
tal mode (determined by the values m, n) is small; see 
Figs. [2] and [3] In this regime the signal and idler photons 
propagate essentially as plane waves, yet with a modified 
dispersion relation with respect to an equivalent source 
without a photonic crystal structure. From Eq.(|7]) there 
is no second-order band gap for a 0.5 duty cycle, which 
implies that the effect of the photonic crystal tends to 
be insignificant for the pump field, thus suppressing con- 
tributions with 1^0. Thus, in this paper, we will con- 
centrate on the contribution I = m = n = 0, which cor- 
responds to the fundamental, forward-propagating mode 
for the three interacting fields. 



Let us note that for type-II PDC in standard nonlinear 
crystals it is typically sufficient to consider a power series 
expansion of the phase mismatch up to the group velocity 
terms. However, for NLPCs with the signal and idler fre- 
quencies in proximity to one of the band gap boundaries, 
GVD and higher-order dispersion terms are strongly en- 
hanced with respect to an equivalent source without a 
photonic crystal structure. Thus, it becomes necessary 
to consider higher-order terms; here we have considered 
up to quartic order in the frequency detunings. 

In order to facilitate an analysis of the conditions under 
which it is possible to generate two-photon states that 
are close to factorable, we write down the joint spectral 
amplitude in terms of Gaussian functions. This can be 
achieved by modeling the pump envelope function as a 
Gaussian function, 



ap{LLis + uji) = exp [-(ws + - 2wo)^/(T^ 



(21) 



IV. CONDITIONS FOR FACTORIZABILITY 

In order to carry out an analysis of the relationship 
between the various experimental parameters and the re- 
sulting spectral entanglement properties, it is convenient 
to express the phase mismatch as a power series in the 
frequency detunings VsA = ^s,i — ^o, where lOq is the 
degenerate frequency and L is the crystal length 



LAK w iAX(°) 
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+2Tf\2v1 + 2,u,u, + 2vl)\v,v, + i?(5).(17) 

Here, 'i?(5) represents fifth, and higher order terms in 
the detunings while A/C*^") represents the frequency- 
independent term, which vanishes when perfect phase 
matching occurs at (jJq, 



(18) 



Equation( 17 1 is written in terms of the mismatch in the 



jth frequency derivative of the wave number between the 
pump and the signal and idler wave packets /, and a 
term proportional to the jth frequency derivative of the 
pump wave number r, 



(19) 



in terms of 



= {l/]\){d'Kjduj^)\^^^ 
K^^^ = {l/jl){d^Kp/du;n\.=2. 



(20) 



and by approximating the sine function in the phase- 
matching function as a Gaussian funciton i.e. sinc(a;) « 
exp(— 70:;^) with 7 « 0.193. Under these approximations, 
the joint spectral amplitude can be expressed as 



f{Vs,Vi) 



exp{-(7/4)[L2AX2 + ^(^^ + 

+iLAK/2}. (22) 



By writing LAK in Eq. ( 22 1 as its power series expan 



sion [Eq.( 17 1] , and keeping terms up to fourth order while 
assuming that Aif'^*'^ — 0, the joint spectral amplitude 
becomes 



/(l/,, J/i) « exp {"(7/4) [$^(i/,) + + $^,(j^s, i^»)]} 

(23) 

where 



4 



2r(iV(^V' - *(2/7) E rk"i^' (24) 



and 



^s^ilys,ly^)/{2ly,l.^) = ri^Vf ^ -|- 4/(7^^) 



-(2rfri^^ 

.(2)^(2) 



3rfV.« 
3^^M^) 



+ [2r^^'T> 

+ [2(rf )2 + rPrf ) + 3rf (ri^^ + rf ))]..., 

-t{2/j)[{4'^ + {3/2)4'\l^s + 'y^) + 
ri")(2z.2 + 3i.,i., + 2i.2)}]. 



(25) 
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While ^fj.ii') (with /i = s,i) represents the contri- 
butions that depend only on one of the two frequen- 
cies, ^sii'^s^'^i) depends on both frequencies, and gives 
rise to correlations in the two-photon state; factorabil- 
ity is attained if ^siit^sj^i) = 0. It is apparent from 
Eqs.(24l and (251 that the dominant terms which govern 



the phase matching behavior are the group velocity mis- 
match terms, proportional to t^V- consider the 
effect of making both of these group velocity mismatch 
terms vanish, i.e.. 



condition, the expression for the joint spectral amplitude 
is thus further simplified. The modulus may now be writ- 
ten as 



\f{iys,Vi)\ ~ exp[-(i^5 + i^i)^/o-^ 

-(7/4)(rp).,^ + rf .,^)^] (30) 

while the argument may be written as 



-(1) -^(1) - 



0. 



(26) 



In this case, for which the pump, signal, and idler prop- 
agate at the same group velocity, the expression for the 
joint spectral amplitude greatly simplifies, and in partic- 
ular many of the mixed terms giving rise to correlations 
are suppressed. The modulus of the joint spectral ampli- 
tude becomes 



\f{u,,t.,)\ « exp[{-;.2/^' -7(rf V^.V4}] 

Xexp[{-^2/^2„^(^(2))2^4/4|] 

+7[(rf )' + rPrfV2]^>,?}, (27) 
while the argument becomes 



arg[/(i^s,i^»)] 



4 

(1/2) 

J=2 



{3)^1 



+ rr^^i\). (31) 



Let us note that while the modulus contains mixed 
terms proportional to Vg^i and to v^vf, the argument 
does not contain mixed terms. Let us now assume, in 



addition to group velocity matching [see Eq.(26l] and 



weak pump dispersion (see Eq.(29l] that the pump is 



broadband. It may be shown from Eq.(30l that the joint 



spectral amplitude exhibits no dependence on the pump 
bandwidth if 



a»2(4/7)i/4(r?^+Tf' )-i/2. 



(32) 



Thus, for a sufficiently broadband pump [so that 
Eq.(32) is satisfied], the modulus of the joint spectral 



amplitude reduces to 



arg[/(i/„ I.,)] « (1/2) E (r^s'^'^i + rf^Vf ) + v,v.\rf^ 



i=2 



+-iT^^\v, + V,) + 2r^^'(2z.2 + Zv^.y^ + 2vl)\. (28) 



[/(:.„ v^)\ « exp[-(7/4)(ri2)^2 ^ 



(33) 



We have seen that when the following three conditions 
are satisfied: (1) complete group velocity matching [see 



Equation (27l is written as aproduct of three exponen- Eq.(26)], (2) weak pump dispersion [see Eq.(29l], and (3) 



tials, where the first two represent the factorable compo- 
nents and the third term describes spectral (modulus) 
correlations between the signal and idler photons. Like- 
wise, in Eq.( |28[ ) it is the second term, with a VsVi overall 
multiplicative factor, which describes spectral phase cor- 
relations between the signal and idler photons. 

Let us now in addition suppose that the signal and idler 
photons experience much stronger dispersion than the 
pump. In particular, let us assume that jth-order pump 
dispersion coefficients Tj, may be neglected with respect 
to the signal and idler jth-order dispersion coefficients 
T^^^ ie 



sufficiently broadband pump [see Eq.( 32 1] , the joint spec 



tral amplitude attains a particularly simple form in which 
there is a single mixed term, proportional to v'^.v'l up to 
fourth order in Afc^. For the specific cases we analyzed 
(see the next section), while conditions 1 and 3 must be 
satisfied in order to attain nearly factorable two-photon 
states, the effect of the mixed terms controlled by con- 
dition 2 is comparatively small. In the next section, we 
will show that the properties of one-dimensional nonlin- 
ear photonic crystals, may be exploited for the fulfilment 
of conditions 1 and 3, and partial fulfilment of condition 
2. 



ri3)\ 



0)1 W^3\ 



(29) 



It may be shown that, if the condition in Eq.(p9| with 



j = 2 is satisfied, the fourth order mixed terms in Eq.(27) 



proportional to may be neglected. Similarly, if the 
condition in Eq. ( 29 1 with j = 2, 3, 4 is satisfied, the mixed 



terms in the argument of the joint amplitude [Eq.(28l] 



may be neglected. Imposing this weak pump dispersion 



V. FULFILMENT OF CONDITIONS WITH 
NONLINEAR PHOTONIC CRYSTALS 

Our strategy for controlling the spectral properties of 
PDC photon pairs generated by one-dimensional NLPCs 
is based on the observation that the group velocity is 
sharply reduced for frequencies in the vicinity of the band 
gap boundaries, as illustrated in Fig.[2](b). It is possible 
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to exploit the properties of photonic crystals to compen- 
sate for material dispersion and to impose specific con- 
ditions on the group velocities of the three interacting 
fields. Indeed, while for a standard optical material the 
pump will tend to propagate at a lower group velocity 
than the generated light, it is possible to design a NLPC 
so that the group velocities of the signal and idler photons 
in proximity to one of the band gap boundaries, are re- 
duced sufficiently to make them equal to the pump group 
velocity. The former can be achieved while maintaining 
essentially a plane-wave character for the corresponding 
signal and idler Bloch waves. Likewise, the NLPC can 
be designed so that the pump frequency is far from band 
gaps, so that the pump field is essentially unaffected by 
the photonic crystal structure. Indeed, NLPCs make 
it possible to attain, even in a type-II geometry where 
each of the three fields experiences different dispersion 
characteristics, complete group velocity matching where 
K'p — K'^ — K[, or in the notation of the previous section 

ri^) = = 0. 

In order to simplify our analysis of realistic NLPCs, we 
first limit the parameter space; we assume that the crys- 
tal is operated at room temperature and assume a duty 
cycle a/ A = 0.5 throughout. Likewise, we regard the 
degenerate PDC frequency lOq as a fixed parameter. It 
is indeed remarkable that, in general, nearly factorable 
PDC photon pair sources based on NLPCs permit the 
specification of an arbitrary central emission wavelength, 
and for that matter an arbitrary material, constrained 
only by the ability to fabricate the crystal with the re- 
quired periodic index of refraction variation. Thus, for 
degenerate coUinear PDC at a given freely specified cen- 
tral frequency, we are left with the following experimental 
parameters: Bragg period (A), permittivity contrast (a), 
crystal propagation angle (Opm), crystal length (L), and 
pump bandwidth (cr). 

In what follows we present a numerical analysis for 
the simultaneous fulfi lme nt of basic phase matching 
[Aif'"^ = 0; see Eq.(18l] and complete group veloc- 

for a one-dimensional NLPC 
Let us note that these condi- 



ity matching [see Eq.(26l 
based on a BBO crystal, 
tions are independent of crystal length and pump band- 
width. Thus, we have three variables: Bragg period 
(A), permittivity contrast (a), and crystal propagation 
angle {9pm) with which to satisfy three conditions: (i) 
AX(o) = (ii) K'p = K'„ and (iii) K'^ = K'^ (where the 
pump frequency derivatives are evaluated at 2iOo while 
signal and idler frequency derivatives are evaluated at 

Wo). 

Figure [4J a) shows in {a, 9pm} space, for three differ- 
ent values of A, the contour defined by the condition 
AK'^^^ = together with the contour defined by the con- 
dition = K'^. Note that both contours shift as the 
Bragg period A is modified. Similarly, Fig. ^h) shows, 
for the same three values of A, the contour defined by the 
condition AK^^'^ — together with the contour defined 
by the condition K'^ = K'^. There exists a specific value 
of the Bragg period A for which the three contours meet 



a single point on {a, 9pm} space, yielding the specific val- 
ues of the three parameters (a, 9pm, and A) which sat- 
isfy simultaneously phase matching and complete group 
velocity matching. Assuming Ao — 27rc/a;o = 850 nm 
(which we stress may be freely specified), these values 
are a = 0.028, 9pm = 41.1°, and A = 274.9 nm, obtained 
numerically from the intersection of the three resulting 
contours [see Fig. |4](c)]. For this choice of parameters. 
Fig. [ijjD) shows a plot of the group velocities for the or- 
dinary and extraordinary waves. In this plot we indicate 
that, for a type-II interaction where the pump is an ex- 
traordinary wave and the signal and idler are extraordi- 
nary and ordinary, respectively, we indeed obtain identi- 
cal group velocities for the three fields. Let us note that 
to our knowledge no other reported technique permits 
complete group velocity matching in a coUinear, type-II 
nonlinear parametric interaction. 
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FIG. 4: (Color online) (a) Graphical representation of the 
conditions (i) AK = and (ii) K'p — K'^, and (b) the con- 
ditions (i) AA" = and (ii) K'^ = K[ on {a, 9pm} space, for 
constant A [we have used the values A = 277.8 nm (shown 
in blue), A = 274.9 nm (red) and A = 271. 9nm (black)]. 
Each interesection point indicates a solution for simultane- 
ous phasematching and group velocity matching between the 
pump and one of the generated photons, (c) Graphical rep- 
resentation of the conditions (i) AK = 0, (ii) K'j, = K'^ y (iii) 
K'p = K'i on {a, 6pm} space, assuming A = 274. 9nm. The in- 
teresection point indicates a solution for simultaneous phase 
matching and complete group velocity matching. Note that, 
in the shaded area, the wave number becomes complex, and 
therefore light does not propagate, (d) Group velocity vs fre- 
quency for the ordinary and extraordinary waves; the dotted 
(black) lines indicate the resulting complete group velocity 
matching. 

According to the conditions derived in the previous 
section, in order to guarantee factorability we require the 
fulfilment of the weak dispersion condition [see Eq.(29l], 
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in addition to complete group velocity matching. In par- 
ticular, if this condition is satisfied for j = 2, then three 
of the four remaining mixed terms in the joint spectral 
amplitude (to fourth order in the detunings) may be sup- 
pressed. As illustrated by Fig. [2] when the signal and 
idler frequencies are in proximity to one of the band gap 
boundaries while the pump experiences essentially only 
the material dispersion, group velocity dispersion and 
higher order terms associated with the generated light 
will be greatly enhanced with respect to corresponding 
pump quantities. Thus, for the specific experimental pa- 
rameters yielding complete group velocity matching (see 
the previous paragraph), the resulting second order dis- 
persion coefficients are ri^'' — 6.12fs^, rp' — 0.87fs^ and 



0.16fs^. The condition It, 



(2), 



(2), 



\tP^\ will be 



more accurately satisfied the closer the degenerate PDC 
frequency is to the band gap boundaries for each polariza- 
tion. However, note that, because the band gap bound- 
aries for the two polarizations are spectrally distinct, the 
weak pump dispersion condition cannot be satisfied to 
the same degree for both photons. Thus, for this spe- 
cific choice of parameters, while the three mixed terms 

(2) 

proportional to Tp are reduced, they are not suppressed 
perfectly. 

The third condition for factorability derived in the 
previous section, apart from complete group velocity 
matching and weak pump dispersion, is that the pump 
should be sufficiently broadband [see Eq.(24)]. Note that 
while the Bragg period, permittivity contrast and crys- 
tal propagation angle must have specific values in or- 
der to guarantee complete group velocity matching, the 
pump bandwidth and crystal length must fulfil a compar- 
atively soft condition given as the inequality in Eq.(32l. 
Figure [5] represents the resulting joint spectral inten- 
sity \ap{u!s + aJi)0ooo('^s, for a pump centered at 
Xp = 2ttc/{2lOo) ~ 425 nm with a full width at half max- 
imum bandwidth of 10 nm, with crystal length L — A 
mm and where the rest of the parameters are as specified 
above, in the context of Fig. [4] Note that in this plot we 
have taken into account the complete dispersion (rather 
than a truncated power series approach) and have not 
used the Gaussian approximation. While these approx- 
imations were essential for our analysis, it is important 
to verify the degree of factorability in a two-photon state 
produced by a realistic source. Figure [5ja) shows the 
pump spectral envelope ap(ws -t- Wi) as given by Eq.(21 1, 
Fig. |5]^b) shows the phasematching function 0ooo('^si "^j) 



as given by Eq.(15l, and Fig. |5jc) shows the joint spec- 
tral intensity. It is apparent from the figure that the 
resulting two-photon state is nearly factorable. Indeed, a 
numerical Schmidt decomposition [23, yields a Schmidt 
number of K fa 1.11 (which could be reduced further by 
moderate spectral filtering). 

As has already been discussed, the effective control 
over group velocity properties permitted by photonic 
crystal structures, which may be exploited to obtain 
nearly factorable two-photon states, enhances higher- 



order dispersion terms. This has the effect that, when 
complete group velocity matching is not fulfilled, the con- 
tours of equal amplitude of the phase matching function 
in LOsjLOi space tend to be highly curved. Therefore, the 
technique presented here is not naturally suited for the 
generation of states exhibiting strict correlation or anti- 
correlation in frequency, unless the bandwidth of interest 
is small enough that the curvature of the phase matching 
function may be neglected. 
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FIG. 5: (Color online) Joint spectral intensity, where disper- 
sive effects to all orders are included and where the Gaussian 
approximation has not been used, (a) Pump spectral enve- 
lope with FWHM bandwidth of 10 nm. (b) Phase matching 
function; the two yellow bands (light gray) indicate the band 
gaps for each of the two polarizations, (c) Joint spectral in- 
tensity, given as the product of the pump envelope and phase 
matching function. 



For practical implementations of this technique, it is 
important to analyze the required tolerances for the ex- 
perimental parameters, in order to ensure factorizability. 
Figure [6] shows a plot of the Schmidt number X as a 
function of the three parameters which must satisfy a 
strict condition in order to attain complete group veloc- 
ity matching. We plot the Schmidt number as a function 
of each of these parameters, while maintaining all others 
equal to their nominal values (which yield the minimum 
value of K). Thus, Fig.[6]ja) shows K vs the Bragg period 
A, Fig.pTb) shows K vs the permittivity contrast a, and 
Fig. l6][c) shows K vs the crystal propagation angle 6pm- 
We define the tolerance in variable x as Ax — X2 — Xi 
where xi^2 (with X2 > xi) are the values where K rises 
to K — \/2Kmin in terms of the value of the Schmidt 
number attained when all variables are equal to their 
nominal values Kmin- Thus, we obtain AA « 1.04 nm, 
Aa « 0.011, and A6'™ « 1.23°. 
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FIG. 6: (Color online) Schmidt number K plotted vs (a) 
Bragg period A, (b) permittivity contrast a, and (c) crystal 
propagation angle 6pm- 
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VI. CONCLUSIONS 

We have studied the generation of photon pairs 
by spontaneous parametric down-conversion in a one- 
dimensional nonhnear photonic crystal, with a broad- 
band pump. We have considered a NLPC based on a 
standard x^"^^ crystal with uniaxial birefringence where 
the dispersion for both polarizations has been altered 
from its natural state so that, throughout the length 
of the resulting crystal, the index of refraction for each 
polarization alternately takes the value consistent with 
material-only dispersion and a value that is larger than 
the former, consistent with a small permittivity contrast. 
We have developed a set of conditions that must be satis- 
fied in order to guarantee a nearly factorable two photon 
state. These conditions are: (i) complete group veloc- 
ity matching, where the pump pulse and the signal and 
idler photons propagate at the same group velocity [see 
Eq.(|26[)], (ii) weak pu mp d ispersion [characterized by co- 



efficients Tp \ see Eq.(19)], relative to signal and idler 



dispersion [characterized by coefficients t^J} , see Eq.(29l], 



and (iii) sufficient pump bandwidth [see Eq.(32)]. We 
have shown that the strongly-modified dispersion in the 
spectral vicinity of the bandgap boundaries in a NLPC 



may be exploited in order to fulfil condition (i), and to 
partially fulfil condition (ii). The essential advantage of 
this technique is that for an arbitrary nonlinear material 
operated at an arbitrary spectral range it is in princi- 
ple possible to design a photonic-crystal structure which 
counteracts the material dispersion so as to permit the 
generation of factorable photon pairs in a type-II degen- 
erate, coUinear geometry which is compatible with wave 
guiding and which permits high efficiency photon pair 
splitting. It is, however, anticipated that the fabrica- 
tion of the appropriate photonic-crystal structure may 
pose the main technical challenge in practical implemen- 
tations. Longer Bragg periods, potentially more eas- 
ily fabricated, could be exploited through higher-order 
bandgaps. It is expected that this work may be impor- 
tant for the development of optimized nonclassical light 
sources for quantum information processing. 
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